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Abstract 

Intermediate-energy heavy-ion collisions can produce a spin polarization of the projectile-like 
species. Spin polarization has been observed for both nucleon removal and nucleon pickup processes. 
Qualitative agreement with measured spin polarization as a function of the momentum of the 
projectile-like fragment is found in a kinematical model that considers conservation of linear and 
angular momentum and assumes peripheral interactions between the fast projectile and target. 
Better quantitative agreement was reached by including more realistic angular distributions, de- 
orientation caused by 7-ray emission, and by correcting for the out-of-plane acceptance. The newly 
introduced corrections were found to apply to both nucleon removal and nucleon pickup processes. 
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INTRODUCTION 



Nuclear spin polarization is a necessary condition for many types of physics experiments 
including studies of nuclear structure, nuclear reactions, fundamental interactions, and con- 
densed matter physics. Many of these experiments require the polarization of radioactive 
ion beams (RIBs). For example, ground state nuclear magnetic dipole (/i) and electric 
quadrupole moments (Q) are critical in the study of nuclear structure, especially as one 
moves away from the valley of stability. Currently, the best reach for such measurements 
has been demonstrated by RIBs polarized in fragmentation reactions. 

Quenched magnetic moment values are observed for the extremely neutron- deficient nuclei 
and ^'^Cu. The small /i value of the proton drip line nucleus (5'p=1.3 MeV) leads to 
a large deviation of the deduced isoscaler expectation value (a) for the ^Li-^C mirror pair 
when compared to observed trends for T = 1/2 nuclei (ill. Shell model calculations were 

n 

shown to better reproduce (a) when isospin mixing was included [2]. Intruder mixing with 
the ground state of caused by proton shell quenching may also account for the deviant 
fi value 3|. ^''Cu represents the heaviest = —1/2 nuclide with a known fi value. The 
magnetic moment is smaller than predicted in shell model calculations that assume ^^Ni is 
a good doubly- magic nucleus [4|. 

The quadrupole moments of the B isotopes also provide evidence for changes in nuclear 
structure. Qi^B) was found to be twice as large as expected from shell model calculations 

n 

[5|]. The larger Q was considered evidence of a proton halo covering a neutron core. Near 
the neutron drip line, the experimental Q moments for odd-mass ^^B, ^^B, and ^''B were 
found to be nearly constant, independent of neutron number Q]. The calculated Q obtained 
with standard values of effective charges overpredict the experimental Q for ^^B and ^''B 
isotopes. The result reveals that both proton and neutron quadrupole polarization charges 
may be reduced near the neutron drip line [3]. 

From these examples, it is apparent that moment measurements play an important role 
in nuclear structure studies. As stated earlier, in order to perform any moment measure- 
ment, the nuclei must be spin polarized. There is a strong dependence of polarization on 
the momentum of the fragment nucleus, therefore, it is crucial to know the magnitude of 
polarization prior to the experiment. While fragmentation reactions are key in producing 
exotic nuclei, these nuclei come with low rates. A figure of merit for measurements involving 
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polarization is P^Y, where P represents polarization and Y is yield. The optimization of po- 
larization with yield is important. Improvements in yield will come with the development of 
new RIB facilities, but while yields are small, the ability to accurately predict the expected 
polarization is required for experimental success. Broader applications for the polarized 
beams will be achieved with new RIB facilities, but at the same time, better predictive 
power is needed. 

In general, the spin orientation of an ensemble of quantum states is usually specified by 
a statistical tensor ^. The tensor can be denoted by Pkq{I), where / is the spin of the state 
and k, the rank of the statistical tensor, takes the values k = 0,1, 2, ...21 if I is an integer, 
and k = 0,1, 2, ...2/ — 1 if / is half-integer, q takes integer values between —k and +k. 
The statistical tensor is related to the distribution of magnetic substates with respect to a 
chosen coordinate frame. It is usually possible to specify the magnetic substate distribution 
by the populations P(m) of the 21 + 1 m-substates, where it is assumed that the P(m) are 
normalized so that Yl -P(^) = 1- 

For an m-state distribution with axial symmetry, only the q = components of Pkq{I) 
are non-zero. Spin polarization is defined in terms of the k = 1 statistical tensor, denoted 
Pio, relative to its value for maximum polarization. Specifically, 

Thus, spin polarization is a measure of the orientation of the total angular momentum 
relative to a fixed axis. 

Spin polarization of projectile-like residues from intermediate-energy heavy-ion reactions 
was first observed at RIKEN [9]. Fragments detected at small angles with respect to the 
normal beam axis were shown to have a polarized spin in the low-intensity wings of the 
momentum distribution. A qualitative description of the polarization mechanism was found 
in a model that considers conservation of linear and angular momentum, and assumes pe- 
ripheral interactions between the fast projectile and target. Fig. [1] presents a schematic of 



and 



so the spin polarization is 



the expected polarization and yield for the nucleon removal process for fragmentation of 
a projectile on a heavy target. A systematic study of the spin polarization following few- 
nucleon removal from light projectiles as a function of energy and target was completed by 
Okuno et al. [l^. This study demonstrated that the relation between the outgoing frag- 
ment momentum and sign of spin polarization depended on the mean deflection angle, O^ef- 
Near-side reactions occur for high-Z targets, where the Coulomb deflection will dominate 
the internuclear potential between projectile and target, giving the polarization dependence 
shown in Fig. [H The nucleon-nucleon potential energy governs removal reactions on low-Z 
targets. Far-side reactions prevail in this case, and the sign of the observed polarization is 
reversed. 

Removal 



>- 





Relative IVIomentum 

FIG. 1: Illustration of nuclear spin polarization produced in a nucleon removal reaction at inter- 
mediate energies, for a high Z target. The yield and polarization curves are given relative to the 
incident projectile momentum. The removal schematic is given in the projectile-like rest frame. 




The spin polarization has a near- zero value at the peak of the fragment yield curve for both 
near- and far-side dominated reactions, since \Odef\ is large. This behavior can be qualitively 
understood from the projectile rest-frame diagram in Fig. [H The removed nucleons have 
momentum K. The 2;-component of the induced angular momentum of the projectile-like 
species is iz = —Xky + Yk^, where X, Y are the localized cartesian coordinates of the 
removed nucleon(s), and k^;, ky are the momentum components of the removed nucleons in 
the reaction plane. If the nucleon removal occurs uniformly in the overlap region, X ~ 
Rq, y ~ 0, then = —Xky. Zero polarization will therefore result when the fragment 
momentum equals the projectile momentum, since ky = in the projectile rest frame under 
these conditions. 
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If nucleoli removal is not uniform in the overlap region, Y and the term Ykx can 
contribute to iz- Such a contribution will only be observed experimentally when \6def\ is 
small. The final scattering angle of the fragment is 6l =Odef+^0, where is the change 
in angle caused by the transverse momentum component of the removed nucleons, AO = 
tan^^(— A;^/p). Here p is the total momentum of the projectile-like fragment. In reactions 
where \Odef\^ 0, it is the transverse momentum component of the removed nucleon(s) that 
"kicks" the fragments to small angles, and the resulting polarization is negative since > 
to give positive and Y < for non-uniform nucleon removal as illustrated in Fig. [H 

A Monte Carlo code was developed [l^ based on the ideas discussed above to simu- 
late the spin polarization generated in nucleon removal reactions at intermediate energies. 
The general behavior of the spin polarization as a function of projectile-like momentum 
was achieved, although a scaling factor of 0.25 was needed to reproduce the magnitude of 
polarization observed experimentally. 

Spin polarization via nucleon pickup reactions at intermediate energies was first demon- 
strated at NSCL III]. Positive spin polarization was determined for ^^K species collected 
at small angles in the reaction of '^^Ar projectiles on a ^Be target at 150 MeV/nucleon. 
Fig. [2] illustrates the spin polarization and yield from nucleon pickup reactions. The key 
to understanding the observed spin polarization in the pickup process is that the picked-up 
nucleon must have an average momentum equal to the Fermi momentum oriented paral- 
lel to the beam direction. Souliotis et al. 12|] reached this conclusion from the observed 
shifts in the centroids of the momentum distributions for one and two-nucleon pickup prod- 
ucts. The average projectile-like momentum < p > was found to satisfy the relation 
< p >=< Pp > + < Pt >, with < pp > the average momentum of the incident particle 
and < Pt > the average momentum of the picked up nucleon equal to the Fermi momentum. 
In the rest frame of the projectile-like species, the momentum of the picked-up nucleon will 
be antiparallel to the incoming projectile momentum. 

The 2;-component of orbital angular momentum induced by the nucleon pickup process is 
£z = RAp, assuming a peripheral interaction where the nucleon is picked up to a localized 
position on the projectile given by R in Fig. [2l The spin polarization will be zero when 
the momentum of the picked-up nucleon matches the momentum of the incoming projectile 
[Ap = 0). The zero crossing occurs at the projectile- like momentum p = [{Ap + 1)/Ap\pp, 
where Ap and Pp are the mass number and momentum of the projectile, respectively. A 
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linear increase in 1^ is expected with a decrease in the momentum of the outgoing pickup 



product. Groh et al. 



ll| found that proton pickup reactions follow the trend shown in Fig. 



m except for the low momentum side of the momentum distribution. At low momentum 
values of the pickup products, the matching conditions for pickup will no longer be satisfied, 
and the spin polarization is observed to rapidly approach zero. 

Turzo et al. showed that neutron pickup reactions at intermediate energies behave in a 
similar manner [l3|. They extended the Monte Carlo simulation of Ref. to include 
nucleon pickup and the momentum considerations discussed by Groh et al. Qualitative 
agreement of the observed spin polarization as a function of the projectile-like product was 
realized, as was the case with nucleon removal reactions. But again, a scaling factor of 0.25 
was needed to reproduce the magnitude of the observed spin polarization. Scaling factors 
of the same magnitude required for both nucleon removal and nucleon pickup suggest that 
the same quantitative correction factors should apply to both. 

The kinematical model proposed by Asahi et al. has been successfully employed to qualita- 
tively explain spin polarization at intermediate energies for both nucleon removal and pickup 
reactions. This paper introduces additional considerations to the kinematical model aimed to 
improve the quantitative agreement of the Monte Carlo simulation with the experimentally- 
observed spin polarization in nucleon removal and pickup reactions. 



Pickup 




I 



H h 





Relative Momentum 




K = (kx,ky,k2) 



FIG. 2: Illustration of nuclear spin polarization produced in a nucleon pickup reaction at inter- 
mediate energies. The yield and polarization curves are given relative to the incident projectile 
momentum. The pickup schematic is given in the projectile-like rest frame. 
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MONTE CARLO SIMULATION 

The simulation code developed by Okuno et al. was used to compute the momenta of 
the removed nucleons in Monte Carlo fashion. The kinematical equations from Ref. ^ were 
then applied to calculate the spin polarization (iz/L) as a function of fragment momentum. 
The nucleon removal positions X, Y on the projectile surface are determined by the relations 
X = Rq cos and Y = —Rq sin 0, where O > as shown in Fig. [H Improvements made to 
the original simulation code, which advance toward better quantitative agreement with spin 
polarization in both nucleon removal and nucleon pickup reactions, are discussed below. 



Location of Nucleon Abrasion 



One shortcoming of the original simulation code was that the projectile-target interaction 
was one-sided. We have adopted an absolute coordinate system to maintain the correct 
relationship between the sign of the polarization, the emission angle of the fragment, and 
the momentum of the fragment. Interactions between projectile and target are permitted in 
the reaction plane, on both "sides" of the target, as shown in Fig. [31 The coordinate system 
is defined as follows. Positive y is the beam direction, positive x is defined to the right of 
the target relative to the beam direction, and positive z is perpendicular to the scattering 
plane forming a right-handed coordinate system. Positive angles are defined to the left of the 
y-axis, or toward negative x. For left-sided interactions, a near-side (far-side) collision will 
scatter to the left(right), or to positive (negative) angles. For the interactions on the right 
side of the target, a near(far)-side interaction must scatter to the defined negative (positive) 
angles, thus the signs of the mean deflection angles are changed for these events. 

One of the fundamental phenomena in fragmentation reactions related to polarization is 
that no frag ment spin polarization is observed when the mean fragment angle, 6l, is zero 
degrees [14]. Shown in Fig. H] is the calculated polarization for 109.6 MeV/nucleon ^^N 
fragmented in a target of ^^''Au to make ^^B at a mean fragment angle of ° with an angular 
acceptance of 0.25 °. The calculated polarization is equivalent to zero within statistical error. 
Allowance for right- or left-sided collisions correctly accounts for the observed absence of 
spin polarization. On average, half of the events detected at ° are from a right-sided 
interaction, while the other half are from left-sided interactions. The resulting polarization 
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FIG. 3: Schematic diagram of right- and left-sided colhsions. The dotted hnes are far-side interac- 
tions and the sohd hnes are near-side interactions. 
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FIG. 4: (Color online) Calculated polarization as a function of relative fragment momentum for 
the fragmentation of 109.6 MeV/nucleon ^^N in a target of ^^''Au to make ^^B at a fragment angle 
of 0°. Po the momentum of the incident projectile and p is the momentum of the outgoing 
fragment. The grey band represents the range of simulation results within a la distribution. 

for each type of interaction is equal in magnitude, but opposite in sign, giving an average of 
zero polarization. 

Previously, an average value for the location X, Y of the removed nucleons in the projectile 
was calculated based on the rotation angle G. The same removal location was then used for 
every event in the simulation. The removal location is now calculated based on the volume 
of intersection of a cylinder and sphere following the prescription of Gosset et al. [15|. The 
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number of nucleons removed is directly proportional to the overlap volume. A removal 
position was calculated in Monte Carlo fashion for each nucleon in the overlap region and 
then averaged to give the position X, Y of the group of removed nucleons. The rotation 
angle is accommodated in this approach by offseting the Y position by a fixed value, yet 
ensuring that the offset Y value remains in the overlap region. 

Distribution of Deflection Angles 

The mean deflection angle ddef is calculated by way of numerical integration. Required 
input parameters are the Z and A of both projectile and target, the energy of the projectile, 
the distance of closest approach, r^m, and the real part of the optical model potential for 
the nucleus-nucleus interaction, Vq. The relations of Gossett et al. [l^ are used to determine 

In the original incarnation of the simulation code, all resulting fragments scatter to a 
single angle defined by Odef-, based on the number of nucleons removed plus the angular 
impulse given to the fragment in the transverse direction. However, for a given number 
of nucleons removed, fragments will scatter to a range of angles whose mean is the mean 
deflection angle. To better simulate experimental conditions, a distribution of scattering 
angles whose average is Odef has been considered, including Rutherford, Fermi, and step 
function distributions. The Rutherford distribution describes classical elastic scattering and 
is strongly peaked at °. This distribution was problematic because fragmentation reactions 
are not classical elastic scattering, and the Rutherford distribution is incompatible with a 
mean angle larger than about 1 ° . The step function distribution ranged in angle from zero 
to twice the mean deflection angle, with the functional form shown below: 



Even though the mean is Odef-, the angular distribution is non-physical. In particular, the 
maximum number of fragments did not occur at 0° as observed in experiment. The Fermi 
distribution (i.e. a Woods-Saxon shape) has the form 






(5) 
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where a = 0.5, 6def is the scattering angle and the value of x is varied to make the mean of 
the distribution equal to the mean deflection angle. The correct mean for the distribution 
could only be achieved with a negative value for x, which is unphysical. 

A compromise was made on a straight-line distribution with negative slope. This distri- 
bution has the advantage of peaking at 0°, and by changing the slope, the desired mean 
deflection angle can be imposed on the distribution. A calculated angular distribution with 
^de/=3 ° is shown in Fig. [51 
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FIG. 5: Linear angular distribution with 0rfe/=3°. 
Out-of-Plane Scattering 

In the original simulation code, projectile-target interactions that occurred only in the 
X, y plane were considered. However, fragments should interact with the target on the top, 
the bottom, and all the angles in between, not just on the right or left side of the target in the 
horizontal x-y plane. Such non-equatorial interactions would decrease the z-component of 
angular momentum in the equatorial (reaction) plane, as shown in Fig. [61 The primed frame 
represents non-equatorial scattering and can be represented as a rotation of the coordinate 
system about the y-axis through the angle (3. 

The resulting ^-component of angular momentum for the non-equatorial scattering is 
L'^ = Lz cos f3 — Lj. sin /5, where is the component of angular momentum projected on the 
X axis and (3 is the rotation angle. In theory, f3 can range from —7t/2 to +7r/2, since the 
projectile interaction can occur anywhere in the yz plane. Integrating all possible contribu- 
tions from the out-of-plane acceptance, and including a l/vr normalization factor from the 
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interval of integration gives 



1 p 2 

= — [Lz cos (3 — Lx sin f5)d(3 = —L^ 

TT / IT 71 



(6) 



Experimental devices have a limited angular acceptance, and thus the range on (3 will 
be smaller than — 7r/2 to +7r/2. Fig. [7] shows a schematic of the beam view of an angular 
acceptance of 2.0°±0.5° horizontally and ±2° vertically. The out-of-plane acceptance of 
±2 ° limits the range on (3 to -53 ° to +53 ° for this case. Equation [6] with these limits 
of integration and a re-normalization of 106 ° rather than 180 ° gives = 0.86L,. As the 
acceptance window moves further from the beam axis, the range on the angle (3 will decrease. 
This causes the value of L'^ to approach that of L^. Various example acceptance windows, 
P angles (i.e. ranges on l3) and the resulting corrections for in terms of L'^ are listed in 
Table [H A vertical acceptance of ±2 ° is assumed in all cases, common in most fragment 
separators. 

A factor of 2/7r has been included in the simulation code as a multiplicative factor on 
the polarization. As shown in Table [U different angular acceptance windows could slightly 
increase this value. 
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FIG. 6: (Color online) Representation of non-equatorial scattering of fragments (primed axes) and 
the scattering of fragments purely in the horizontal x — y plane (unprimed axes). Labels defined 
in the text. 
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7-ray Deorientation 



The fragmentation process will often leave the projectile-like product in an excited state. 
One pathway to remove the excitation energy is through 7-ray emission. Discrete 7 rays have 
been observed in "in-beam" fragmentation studies [l6j. Intermediate-energy reactions also 
populate high-spin isomeric states, whose depopulation can be readily followed in the absence 



of the prompt radiation background 17|. 7 rays emitted from a spin polarized nucleus can 
reduce the magnitude of the polarization. Such a deorientation process was not previously 
considered. 7-ray deorientation was included in the simulation code by assuming a statistical 
cascade through a continuum of levels following a prescription similar to that of Leander 




FIG. 7: Representation of the angular acceptance window from a beam view. The fragment angle 
acceptance (horizontal acceptance) is 2.0 "±0.5° and the vertical acceptance is ±2°. 



TABLE I: Polarization correction factors due to non-equatorial scattering. 



Fragment angle(°) 


Acceptance( °) 


Pmax ( ) 


Correction factor 


1.0 


±2.5 


90 


0.64 


2.0 


±0.5 


53 


0.86 


2.0 


±2.5 


90 


0.64 


4.0 


±2.5 


53 


0.86 


5.0 


±2.5 


39 


0.92 





full 


90 


0.64 
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18 1 . A nuclear level density is specified by a constant temperature level-density formula 
191]. For each 7 ray, the transition energy, multipolarity {E1,E2, or Ml), and spin change 
are determined by random numbers. The deorientation coefficients Uk [20I] are calculated 
from the initial and final spin values, and the multipolarity. A cumulative deorientation is 
obtained from the product of the Ui values of the cascade 7 rays. The statistical nature of 
the cascade decay is included by repeating the calculation with a random walk. The resulting 
average deorientation coefficients are given a weighting factor based on the calculated spin 
distribution and excitation energy of the fragment. 

Deorientation coefficients were calculated for the fragmentation reactions involving two- 



nucleon removal 
Refs 



reported in Ref. [lO|, as well as the nucleon pickup reactions reported in 



12] and 131]. The values are presented in Table [TTl The Ui values ranged from 
0.50 to 0.89, suggesting that 7-ray deorientation can have significant impact on the final 
spin polarization of the projectile-like fragment. The calculated average entry spin (~ 2/i) 
and cascade multiplicity (~ 3) were reasonable based on in-beam spectroscopy of light 
fragmentation products |l6[ ]. 



TABLE II: Polarization correction factors due to 7-ray deorientation 



Reaction 


Energy (MeV/A) 


Ui 


Avg. Entry Spin 


Multiplicity 


197Au(14n,12b)X 


39.4 


0.500 


1.54 


2.32 


19Uu(15N,13b)X 


68 


0.592 


1.79 


2.35 


197Au(15n,13B)X 


109.6 


0.576 


1.71 


2.34 


93Nb(i5N,i3B)x 


67.3 


0.583 


1.74 


2.34 


27a1(15n,13b)X 


68 


0.555 


2.02 


2.39 


9Be(36Ar,37K)X 


150.0 


0.688 


2.59 


2.95 


9Be(36S,34Al)X 


77.6 


0.887 


4.17 


2.82 



It can be anticipated that nuclear structure effects, especially the discrete level sequence 
at low energy, may affect the magnitude of deorientation in some cases. For example, ^^Be 
has an = 1/2^ ground-state and its only excited state at 320 keV has I'^ = l/2~. Any 
population of the ground-state via the excited state will have the opposite polarization to 
the direct population of the ground-state. This is an exceptional case. In p-shell nuclei 
closer to the valley of stability that have low lying 1/2+ and 1/2" states, with one being 
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the ground state, the excited spin 1/2 state is accompanied by several other states of higher 
spin that decay directly to the ground state. Thus, if the population of the excited states 
is not selective, the effect of the change in the sign of the polarization for decays through 
excited spin 1/2 states is not expected to be prominent in most cases. 

Low-lying spin states also have the potential to greatly diminsh the net polarization. 
For example, in ^^N, which has a 2~ ground state, the first-excited state at 120 keV is a 0~ 
state. Any population that passes through this state will lose its orientation. However, in 
^^N there are 3~ and 1~ states at 298 and 397 keV, respectively, which will compete with the 
0~ state for population as a result of the reaction, and ameliorate the overall polarization 
loss. 

Thus although specific nuclear level sequences at low spin are important, especially in 
weakly bound systems with very few excited states, as a first estimate of the magnitude of 
the effect it seems appropriate to model the deorientation in a generic way by assuming a 
statistical cascade through a continuum of levels using a Monte-Carlo simulation as described 
above. 



Nucleon Pickup 

The simulation code has been modified to include nucleon pickup, independent of the 



efforts reported in Ref. ISj. The pickup process follows the observations of Souliotis et 



al. [12|], in that the picked- up nucleon has an average momentum equal to the Fermi mo- 
mentum (230 MeV/c), oriented parallel to the beam direction. The momentum distribution 
for the one-neutron pickup reaction ^''A1(^^0,^^0) at 80 MeV/nucleon shows a considerable 
shift of the centroid below the momentum/nucleon of the beam, as observed in Ref. 



3, 



m 



contrast to those for nucleon removal products. The simulated position of the centroid agrees 



with the calculation of Ref. 121], where a simple model based on momentum conservation 
was used (see Fig. [8]). 

The width of the momentum distribution is observed experimentally to be small (around 



20 MeV/c), while it is calculated to be zero, based on the equation by Goldhaber [2l|, and 
extended to nucleon pickup: 

^11=^0 ^^31 (7) 

where App = Ap — AA^ is the mass of the projectile part of the final product and AAt is the 
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FIG. 8: (Color online) Parallel momentum/nucleon distribution calculated with the simulation 
code for the reaction ^^Al(^^0/^0) at 80 MeV/nucleon. The red squares are the data points and 
the black line represents the simulation results. The arrow corresponds to the momentum/nucleon 
of the beam. The simulated momentum distribution has been scaled by the ratio observed in Ref. 
[l^ of experimental centroid/calculated centroid (0.969/0.978), in order to compare to the data. 

number of nucleons picked up from the target. The parameter cxo is the reduced width, and is 
related to the Fermi momentum of the nucleon motion inside the projectile: cTq = VFermJ^- 
A value of ctq = 80 MeV/c is used, which agrees with the experimental widths in Ref. 12|. 
Eq. [7] assumes that the nucleon is picked up from the target with a fixed momentum and 
direction, and the picked up nucleon makes no contribution to the width. Thus, for any 
pure nucleon pickup process, Ap = App and the parallel width is zero. In order to model 
experimental observations of Ref. 12[|, a parallel width of 20 MeV/ c was used. In addition to 



the parallel width. Van Bibber et al. 



22j showed that in heavy-fragment studies in the 100 



MeV/nucleon region, the projectile is subject to an orbital deflection due to its interaction 
with the target nucleus before fragmentation takes place. The orbital deflection gives an 
additional dispersion of the transverse momentum, as demonstrated in the expression: 

o 2^pf{Ap — App) ^App{App — l 

- + CT2- 



(J 1 



a, 



(8) 



Ap-1 ' Ap{Ap - 1) 

where the first term was defined previously, and the second term contains the variance 
of the transverse momentum of the projectile at the time of fragmentation (200 MeV/c as 
used in Ref. {22]). A comparison of the simulated momentum distribution is shown in Fig. 
[8] with the data taken from Ref. [l^ . 
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FIG. 9: Variable definitions for mean deflection angle calculation. 
Optical Potential 

The real part of the optical model potential for the nucleus-nucleus interaction, Vq, is a 
required input parameter for the mean deflection angle calculation. For a single interaction, 
the deflection angle 9 (see Fig. [9]) is given by 

9 = n-24> (9) 

with 

(10) 



where b is the impact parameter, r is the distance between the centers of the two objects, 
U{r) is the potential governing the interaction of the two objects, r^m is the separation 
between the centers of the two point-like objects at the distance of closest approach and the 
energy, E, is given by 

E = ^mvl (11) 



23| 



where Voo is the velocity of the projectile at r = cxo 

The projectile is assumed to move away from the target after the scattering event with 
momentum equal to the incident momentum, thus E{y oo) = E{y incident)- This formula is 
general for any spherically symmetric potential. 

The potential U (r) is defined by 

Uir) = Ucoulomb{r) + Unuclear{r). (12) 

The Coulomb part of the potential is repulsive and is equal to 

1.438Z,Z, 

Ucoulomb{r) = ^ — (13) 

16 



where and Zt are the charges on the projectile and target respectively, and r is the 
separation in fm. The nuclear part of the potential is based on the real part of the optical 
model [2J], and is attractive: 



Here Vq is the depth of the optical model potential, R = 1.2{^/A^ + v^A) where Ap and At 
are the masses of the projectile and target respectively, and a is a measure of the diffuseness 
of the nuclear surface. Vg and a are parameters fit to experimental data. There are very 
limited nucleus-nucleus scattering data, and an exact determination or parametrization of 
Vo is difficult for any given projectile-target combination. Typically this is not a problem 
because in head-on collisions, the nuclear potential does not have a large influence. However, 
the treatment of peripheral collisions depends on the optical potential. In the minimum, a 
determination of Vg is needed. A parametrization of Vq based on energy and/or number of 
nucleons removed would suffice, but unfortunately, such a parametrization does not presently 
exist. 

In the reactions studied in Ref. 



lOj, Vq was calculated by the authors, and the corre- 
sponding 9def was reported in the literature. The reported Odef (given in Table IIIip was used 
in the polarization calculations to follow. 

For the nucleon pickup reactions, Vg was determined with a folding model calculation 



251]. The model was chosen because it reproduces experimental scattering data for heavy 
ions in the energy range of interest. The folding calculation yields the real part of the optical 
potential ( Vq) as a function of the internuclear radius, the distance between the center of the 
projectile and target. The internuclear radius is calculated in the simulation code, based on 



the relations by Gosset et al. [15|, as mentioned previously. The value of Vg corresponding 
to a radius for one nucleon overlap was used in the nucleon pickup reactions. Given in Table 
HIT] are the Odef and Vq for the reactions studied in the present work. 

Results 

Comparison of the spin polarization measurements for intermediate-energy reactions that 
include two-nucleon removal in Ref. lOj with the simulation results that include consider- 



ations for angular distributions, out-of-plane scattering, and 7-ray deorientation, are pre- 
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sented in Fig. [TOl No scaling factor was used to adjust the simulated spin polarization. 

The predicted polarization trend as a function of outgoing fragment momentum for a 
three-nucleon removal case is shown in Figure [TT] The polarized ^^Ni nucleus from a ^^Ni 
primary beam is proposed to determine the ground state magnetic moment measurement 
of ^^Ni at NSCL |28|. The expected rate of ^^Ni is 500 pps, thus optimized polarization at 
peak yield is critical to realize a successful measurement. The proposed measurement will 
complete the fi measurements of nuclei one nucleon removed from ^^Ni. The two Tj = +1/2 
nuclei, ^^Co and ^''Ni, both have known magnetic moment values that compare favorably 
with results of shell model calculations. On the other side of N = Z, the = —1/2 nucleus 
^'^Cu has a /i value that deviates significantly from shell model expectations, suggesting a 
weakened core in ^^Ni. The measurement of /i(^^Ni) is key to understanding how the ^^Ni 
core may be polarized when coupled to a neutron hole. 

The simulation results for the one nucleon pickup processes are shown in Fig. [12] and [T3l 
Angular distributions, out-of-plane scattering, and 7-ray deorientation were implemented in 
the nucleon pickup process as well. In the momentum distribution calculations discussed 

n 

previously, Souliotis et al. used the "typical" Fermi momentum of 230 MeV/c. In order 
for the simulation to encompass a range of targets, the Fermi momentum was taken not 
as 230 MeV/c, but was calculated based on data taken from Moniz et al. 26|. The Fermi 
momentum ranges from 170 MeV/c for the lightest targets to 260 MeV/c for heavier targets. 
The results of the simulation for proton pickup are shown in Fig. [12], calculated for the one 
proton pickup reaction, ^Be(^^Ar,^^K)X, first observed by Groh et al. ll(|. The momentum 



TABLE III: Odef and corresponding Vq used for the reactions studied in the present work. 



Reaction 


Energy (MeV/A) 




Vo (MeV) 


19Uu(14N,12b)X 


39.4 


3.09 


48 


197Au(15n,13B)X 


68.0 


0.92 


63 


197Au(15n,13b)X 


109.6 


0.09 


85 


93Nb(i5N,i3B)X 


67.3 


-0.35 


65 


27a1(15n,13b)X 


68.0 


-3.25 


118 


9Be(36Ar,37K)X 


150.0 


-0.07 


29 


9Be(36s,34Al)X 


77.6 


-0.49 


32 
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(P-Po)/Po 



FIG. 10: (Color online) Polarization as a function of fragment momentum p relative to the 
primary beam momentum po for the removal reactions ^^'^Au(^^N/^B) at 39.4 MeV/nucleon, 
197Au(15n,13b) at 68 MeV/nucleon, 19'^Au(15N/3b) at 109.6 MeV/nucleon, 93Nb(i5N,i3B) at 67.3 
MeV/nucleon, and 27a1(15n,13b) at 68 MeV/nucleon. The red squares are the experimental data 



lot]. The black 



points and the blue dashed lines are the previous simulation results, both from Ref. 
band represents the range of the present simulation results within a la distribution. Momentum 
is given relative to the peak of the yield distribution. 



matching conditions 27| for simple surface-to-surface pickup are best met for the two data 
points on the high momentum side of the yield distribution, where the simulation agrees 
with the data. On the low momentum side of the peak of the yield curve, the picked-up 
nucleon has a momentum less than the Fermi momentum, and the momentum matching 
conditions for direct pickup are poorly satisfied. More complex transfer mechanisms are 
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(p-Po)/Po (p-p„)/p. 



FIG. 11: Polarization distributions generated from the Monte Carlo simulation code as a function of 
relative fragment momentum for the three-neutron removal fragmentation reaction with a primary 
^^Ni beam at 140 MeV/nucleon on a a) ^Be target and b) ^^'^W target. Both targets will be used 
during the experiment, to provide better comparison to the simulation. 



therefore required to describe the polarization on the low momentum side 




{P-PoXPo 



FIG. 12: (Color online) Polarization as a function of fragment momentum p relative to the primary 
beam momentum po for the one-proton pickup reaction ^Be(^^Ar,^^K) (150 MeV/nucleon). The 



111 ] and the grey band represents the range 



red squares are the experimental data points from Ref. 
of simulation results within a la distribution. 

The simulation code was also used to model the neutron pickup data obtained in Ref. [3] , 
as shown in Fig. [131 Again, no scaling factor was required to achieve a better quantitative 
agreement with experiment. 
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FIG. 13: (Color online) Polarization as a function of fragment momentum p relative to the mo- 
mentum at the peak of the yield distribution po for the one-neutron pickup reaction ^Be(^^S,^'^Al) 
(77.5 MeV/nucleon). The red squares are the exjaerimental data points and the blue dashed line 
is the previous simulation result, both from Ref. p^]. The grey band represents the range of the 
present simulation results within a la distribution. 

SUMMARY 



A statistical Monte Carlo code for nuclear spin polarization has been modified to include 
corrections to the kinematical model that aim to improve quantitative agreement with ex- 
periment. The positions of removed nucleons are assigned in a Monte Carlo fashion, and 
projectiles are allowed to interact on either side of the target. The real part of the optical 
potential, Vq, is obtained with a folding model calculation and then used in the calculation 
of the mean deflection angle, Odef- Fragments are permitted to scatter to a distribution of 
deflection angles rather than a single mean angle, and the out of reaction plane acceptance 
has been taken into account. The process of de-orientation due to 7-ray relaxation is also 
included. The angular distribution implementation reduces the calculated polarization mag- 
nitude by about 10% and the out-of-plane acceptance reduces the polarization magnitude 
by about 40%. The calculated polarization magnitude is corrected by about 50% due to 
7-ray de-orientation. These corrections account for the 0.25 scaling factor needed by both 



Ref. 



10[ and Ref. 13j in their simulations. 



It should be noted that the corrections are approximations. A primitive function is used 
to calculate the distribution of deflections angles, and no attempt was made to rigorously 
reproduce the scattering angles observed in fragmentation reactions. The out-of-plane scat- 
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tering correction is implemented as a constant multiplicative factor based on the expected 
acceptances of typical fragment separators. The deorientation correction due to 7-ray relax- 
ation is calculated in a statistical manner considering average excitation energy and entry 
spin. Even within these limitations, the magnitude with which the corrections are shown 
to improve the results indicate that these are important contributions to consider in the 
simulation of polarization in intermediate-energy reactions. 
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